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Abstract 

This paper considers plasmas in which the electrons and ions may have different temperatures. 
This is a case that must be examined because nuclear fusion processes, such as those that appear 
in ICF capsules, have ions whose temperature runs away from the electron temperature. A fast 
charged particle traversing a plasma loses its energy to both the electrons and the ions in the 
plasma. We compute the energy partition, the fractions E(,/Eq and Ei/Eq of the initial energy 
£"0 of this 'impurity particle' that are deposited into the electrons and ions when it has slowed 
down into a "schizophrenic" final ensemble of slowed particles that has neither the electron nor 
the ion temperature. This is not a simple Maxwell-Boltzmann distribution since the background 
particles are not in thermal equilibrium. We perform our calculations using a well-defined Fokker- 
Planck equation for the phase space distribution of the charged impurity particles in a weakly to 
moderately coupled plasma. The Fokker-Planck equation holds to first sub-leading order in the 
dimensionless plasma coupling constant, which translates to computing to order nlnn (leading) 
and n (sub-leading) in the plasma density n. An examination of the energy partition for the 
general case, in which the background plasma contains two different species of particles that are 
not in thermal equilibrium, has not been previously presented in the literature. We have new 
results for this case. The energy partitions for a background plasma in thermal equilibrium have 
been previously computed, but the order n terms have not been calculated, only estimated. Since 
the charged particle does not come to rest, but rather comes into a statistical distribution, the 
energy loss obtained by a simple integration of a dE/dx has an ambiguity on the order of the 
plasma temperature. Our Fokker-Planck formulation provides an unambiguous, precise definition 
of the energy fractions. For equal electron and ion temperatures, we find that our precise results 
agree well with a fit obtained by Fraley, Linnebur, Mason, and Morse. The "schizophrenic" final 
ensemble of slowed particles gives a new mechanism to bring the electron and ion temperatures 
together. The rate at which this new mechanism brings the electrons and ions in the plasma into 
thermal equilibrium will be computed. 
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I. INTRODUCTION 



The underlying theme of this paper is the thermonuclear burn of deuterium-tritium plas- 
mas. We do not consider the initiation of the burn process, which is system specific, nor are 
we interested in the late stages of the process when most of the DT fuel has been burned into 
alpha particles and neutrons, and the electrons and ions arc nearly in thermal equilibrium. 
We instead focus on intermediate times when, in general, there is a significant difference be- 
tween the electron and ion temperatures, but the alpha particle density has not yet become a 
significant fraction of the D and T ion densities.^ The fusion rate is very sensitive to the ion 
temperature Tj. The ion temperature is determined by competition between deposition of 
the alpha particle energy into the ions, which of course increases Tj, and thermal equilibra- 
tion with the electron distribution, which drives Ti down. Our main concern in this paper is 
the partition of the total alpha energy between the ions and electrons in a two-temperature 
plasma in the circumstances that we have outlined.^ This is important in the understanding 
of the time scale and the robustness of the fusion process. Our evaluations of the functions 
which determine the energy partition do not include a contribution from the alpha particles; 
hence our results are valid only if the ensemble of alphas is sufficiently dilute. We find that 
the alpha particles slow down into a non-Maxwellian distribution in which the mean alpha 
energy E lies between the thermal energies of the ions and electrons. Our work shows that 
these non-thermal alpha particles increase the rate of energy transfer between the electrons 
and ions but, since we do not examine late times where the population of alpha particles 
is large, this new mechanism docs not significantly enhance the energy transfer rate. In 
general, as in other work on stopping power and the partition of a fast impurity particle's 
energy to the electrons and ions in the plasma, we assume (as is most often the case) that 
the stopping times are much shorter than the time scale of the fusion so that we can work 
in the adiabatic approximation in which the time dependences of our results are only those 
brought about by the changes in the plasma parameters on which they depend. We also 
require, as is also generally assumed, that the charged particle range is short in comparison 
with the distances over which the plasma conditions vary so that the plasma may be treated 
as being uniform. 

The major results of this paper are as follows. First, as we have mentioned in the previous 
paragraph, we have worked out the energy partition for differing electron and ion tempera- 
tures; this has not been previously considered in the literature. Second, even for the case of 

^ When the alpha particle density is a significant fraction of the plasma ion density, the effect of the 
alphas on the dielectric response of the plasma must be taken into account. This introduces additional 
complications, and as such merits a separate publication. 

^ A short preliminary account of the methods that we employ in this paper, but restricted to the case of 
equal ion and electron temperatures, has previously been presented in [1]. 
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equal ion and electron temperatures, where the alphas relax into a Maxwellian distribution, 
we have made two improvements. We have developed a formulation that precisely defines 
the energy partition so that a correction of order T/Eq is now included, a correction that 
is missing in the literature. In addition to the well-known nlan {n is the number density) 
terms in the energy partition, we have computed exactly the coefficient of the order n term, 
which has previously been only estimated. We turn now to describe our work in some detail. 

When a fast charged particle with initial energy Eq traverses a plasma, it loses its energy 
at a rate dE/dx per unit of distance, and it comes into a quasi-static equilibrium state 
after depositing its initial energy into the electrons and ions that make up the plasma. In 
the thermonuclear fusion process of deuterium and tritium, D -\-T ^ n -\- a, which occurs 
in inertial confinement fusion experiments, the amount of the initial alpha-particle energy 
Eq — 3.54 MeV that is transferred to the D, T ions is crucial because a high ion temperature 
is necessary for the fusion reaction parameter {(tv)^ to become sufficiently large so as to 
have a robust and stable fusion burn. 

In the picture in which the projectile traverses linearly through the plasma until coming 
to a complete stop, the energy partition into ions and electrons is given by 

and 



'o ^ Jo dE/dx 
Here dE^/dx and dE^/dx are the stopping power contributions from the ions and electrons, 
and dE/dx is the total stopping power, 

dE _ dE, ^ dEe ,^ 
dx dx dx ' 

and thus 

E, + Ee = Eo. (1.4) 

This simple picture, however, is only an approximation. For a plasma with equal ion and 
electron temperatures, a fast charged particle does not simply come to rest in the plasma, but 
rather, it becomes thermalized at the ambient plasma temperature T — Ti — Te. Expressing 
temperature in energy units, as we shall do throughout this paper, the correct electron-ion 
energy partition relation should read 

E, + E, + ^T = Eo. (1.5) 

Consequently, rather than extending the lower limits of the integrals (1.1) and (1.2) down 
to zero energy, lower limits of order the temperature, S'min ~ T, must be chosen. The two 
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integrals (1.1) and (1.2) have somewhat different thermal cutoffs, both of order T, and this 
simple picture has a systematic error of relative order T/Eq. We see that the correction 
becomes more important as the plasma temperature is elevated. To account for the energy 
partition in a precise fashion, we shall employ the Fokker-Planck equation in the version 
introduced by Brown, Preston, and Singleton (BPS) [2]. We shall find that the correct 
expression for the energy partition does not, in fact, involve the stopping powers dE^/dx 
and dEe/dx, but rather certain ion and electron functions Ai and Ae that enter into this 
Fokker-Planck equation. In the notation of BPS, the stopping power of a particle of energy 
E — ^mv^ is oi the generic form 



dEg 
dx 



^ T d 
mv dv 



:i.6) 



The functions Aj and Ae thus approach dEj/dx and dE^/dx at high energies, but differ from 
these stopping powers at low energies on the order of the thermal background temperature. 

For the case of equal electron and ion temperatures, the explicit evaluations for E^ and E^ 
derived in Eqs. (4.38) and (4.39), which omit of negligibly small exponential terms involving 
exp{— /3£^o}) read: 



El 



and 



Ee = 



A(E) 



A(E) 
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where 



A{E) = A,{E) + Ae{E) , (1.9) 

and (5 = 1/T. Here erf(x) is the error function defined in Eq. (4.37). Using the definition 
(4.37), partial integration can then be used to show that the sum rule (1.5) follows from 
Eqs. (1.7) and (1.8). Since dEh/dx — )■ Ab for large energies, and since the error function 
approaches unity at large /3E, the precise results (1.7) and (1.8) approach the more intuitive 
but less accurate forms (1.1) and (1.2). Significant differences occur only for ~ T. 

We have numerically evaluated these integrals using the expressions for the A functions 
derived in BPS that are reproduced in Appendix A. We shall compare our results with 
the less precise but well known results of Fraley, Linnebur, Mason, and Morse (FLMM) [3] . 
Starting with a model of the stopping power in an equimolar DT plasma, these authors show 
that the simple rule 

1 
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Eo 1 + Tc/Te 
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provides a good fit to their calculations. The crossover temperature Tc, where the electron 
and ion fractions are equal, can be determined from their Fig. lb. Fraley et al. find 
Tc = 32keV at the density p = 0.213 g/cm^, or a corresponding electron number density 
He = 5.0 X 10^^ cm~^. At the number densities rie = 1.0 x 10^"^ , 1.0 x 10^^ , 1.0 x 10^^ cm~^, 
we find, by fitting our more precise results, that Tc = 31 , 30 , 28 keV respectively. Figure 1 
shows our result (1.7) for the fractional energy loss to ions and the FLMM fit (1.10) for a 
DT plasma with electron number density rie = lO^^cm"^. In this comparison, we use the 
more accurate value Tc = 28 keV. In Fig. 2 we compare the differences between our result 
(1.7) and the FLMM fit (1.10) over a wide range of densities. We see that the FLMM fit 
somewhat overestimates the energy deposited to ions for temperatures above 120 keV over 
a wide range of densities. 




200 



T (keV) 

FIG. 1: The fractional energy loss into ions as a function of the plasma temperature for an a particle 
in an equimolar DT plasma with initial energy Eq = 3.54 MeV. The electrons and ions have a common 
temperature T and the electron number density of the plasma is Ue = 1-0 x 10^^ cm~^. The solid red line is 
the evaluation of Eq. (1.7) while the dashed blue hue is the FLMM fit (1.10) with Tc = 28kcV rather than 
the 32 keV value used by FLMM. 

As Figs. 1 and 2 show, the FLMM fit (1.10), modified slightly to use better values of the 
crossover temperature Tc, is in good agreement with our precise results in the case of equal 
temperatures so long as these temperatures are less that about 120 keV. However, as Fig. 3 
demonstrates, this simple form fails to provide an accurate estimate of the energy partition 
when the ion and electron temperatures are significantly different. These results for differing 
electron and ion temperatures follow from Eqs. (4.48) and (4.50). They are spelled out in 
more detail in the tables presented in our concluding section V. 

Although the work of FLMM continues to be used, a more recent evaluation of the energy 
partition has been carried out by Li and Petrasso [4]. Comparing their Table 1 (Ti = T^) 
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FIG. 2: Differences between the fractional energy losses Ei/Eq as given by the precise result (1.7) and the 
FLMM fit (1.10) for an alpha particle with initial energy Eq = 3.54 MeV in an equimolar DT plasma. The 
different curves correspond to the three electron densities rig of 1.0 x 10^^ cm~^ (red), 1.0 x 10^^ cm""^ (blue), 
and 1.0 X lO^^cm"'^ (green), with the fit (1.10) evaluated with Tc = 31 , 30 , and 28keV, respectively. 

with our Fig. 1 shows that their results for Ei/Eq are ~ 10% too high.^ This discrepancy is 
of order the sub-leading corrections to the Coulomb logarithm.^ 

The emphasis in this paper is the energy partition for unequal electron and ion tem- 
peratures, which is of the utmost importance for DT burn since there the ion temperature 
generally runs away from the electron temperature once the fusion process begins. To de- 
scribe this in a simple fashion, we assume that the sources of the charged impurity particles 
(the a particles in DT fusion) are uniformly distributed throughout the plasma, and that 
the particles are emitted isotropically; hence, the phase space distribution of the impurity 
particles is only a function of the energy and time. The evolution of this distribution is gov- 
erned by a Fokker-Planck equation that involves the coefficient functions Ai{E) and Ae{E) 
which were computed in BPS to order n{\nn + c) in the plasma density n. Since n ~ g'^, 
with g the plasma coupling constant, it is evident that these two terms in the density are 
the leading and first sub-leading terms in the perturbative expansions in g of the coefficient 
functions. Higher-order terms in the expansions become significant at high densities, hence 
our results are not apphcable, in particular, to (strongly coupled) warm dense plasmas. Nu- 
merical simulations provide the only potentially reliable means of validating our analytic 

^ A detailed discussion of the results of Li and Petrasso [4] for the stopping power dE/dx was presented in 
the BPS paper [2] that provides the basis for the work which we perform here. 

^ Long and Tahir [5] have also presented results for the energy partition, but they only compute the separate 
electronic and ionic contributions to the stopping powers, dE^/dx and dEi/dx, as a function of the range 
X for equal temperature background plasmas. They do not present the total energies deposited to the 
electrons and ions, and they also do not present a precise Coulomb logarithm. 
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FIG. 3: The red curves show the energy fractions deposited to the ions as functions of the ion 
temperature Tj for various values of the electron temperature Tg for an a particle with initial 
energy £"0 = 3.54 MeV in an equimolar DT plasma with an electron density = 1.0 x 10^^ cm~^. 
The blue horizontal lines are the energy fractions determined by the FLMM fit (1.10) with our 
value Tq = 30keV. The red curves describing larger values of E^/Eq correspond to increasing 
electron temperatures of Tg = 10 , 30 , 50 , 100 , 200 keV. If the fit (1.10) were exact, the red curves 
would cross the blue horizontal lines when Tj = Tg. This condition is fairly well met except at 
the highest electron temperature Tg = 200 keV where the value of Ei/Eq given by the red curve 
at Tj = 200 keV is considerably smaller than the blue horizontal line. This is consistent with the 
discrepancy at these temperatures shown in Fig. 1. 

expressions for the energy partition in weakly coupled plasmas and evaluating the partition 
in moderately to strongly coupled plasmas, though such computations have not been per- 
formed. Careful, large statistics, molecular dynamics (MD) simulations have been carried 
out by Dimonte and Daligault [6] to investigate electron-ion temperature relaxation over a 
wide range of plasma parameters that span weak to strong coupling. Their MD results for 
the Coulomb logarithm for this process agree with those of BPS [2] for g < 0.2 to within 
the statistical uncertainty of ±5% in the simulations. This indicates the range of validity of 
the Fokker-Planck equation that we use to compute a different result, the energy partition. 

Following a detailed discussion of the Fokker-Planck equation in Sections II A and II B, 
the late-time distribution foo{E) of a 6{t) source of impurity particles, which is needed to 
obtain the electron-ion energy split, is derived in Section II C. In Section III a source is 
slowly turned on and eventually emits particles at a constant rate. The solution f{E,t) 
of the now inhomogeneous Fokker-Planck equation is shown to be the sum of two terms: 
f{E,t) = n(t) foo{E) + f{E), where n(t) is the number density of impurity particles that 
have come into the equilibrium state described by foo{E), and f{E), which describes the 
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transfer of energy to the electrons and ions. The energy losses Eg and Ej to the electrons 
and ions are expressed as single integrals involving the function f{E) [which depends upon 
the ^-coefficients] and the ^-coefficients themselves. The late-time ensemble of impurity 
particles with energy distribution /oo(-E') is not in thermal equilibrium with the background 
plasma, i.e. foo{E) is not a Maxwell- Boltzmann distribution. This ensemble increases the 
rate of ion-electron thermal equilibration above that of the impurity-free plasma. In Section 
IV A we carry out the exphcit construction of f{E). We show how our general results 
for the deposited energy fractions E^ and E^ in the equal temperature case reduce to the 
usual expressions involving dE^/dx and dEg/dx in Section IV B and then describe how these 
approximate results are corrected with our precise formulation in Section IV C. In Section 
IV D we compute E^/Eq and E^/Eq for the general case of different plasma electron and ion 
temperatures in terms of integrals over Ai and Ae- The conclusion V provides a summary of 
our major results including a table of the energy fractions E^/Eq and Eg/Eo for a wide range 
of plasma parameters. At this point, we have finished a logically complete exposition of our 
methodology and results, which is essentially self-contained. However, for those interested in 
supporting details and who may wish to work out the intermediate steps in our calculations, 
we include these details in the Appendices. We provide a review of the A functions that were 
computed in BPS [2] which are needed for the present work in Appendix A, a host of details 
on these functions that include their approximate forms in various regions in Appendix B, 
and an accurate approximation for one of the two multiple integrals appearing in our final 
expressions for Ei and Eg is provided in Appendix C. 



II. FORMULATION OF THE PROBLEM 



A. The Fokker-Planck Equation to Leading and Next-to-Leading Order 

We consider a plasma containing a dilute population of "impurity" particles with a phase 
space density f{r,p,t). For example, in a deuterium-tritium (DT) plasma, the impurities 
could consist of the charged a particles produced from the DT fusion. The problem we shall 
address is the manner by which such impurities reach a quasi-static equilibrium distribution. 
During this process, the impurities deposit portions of their energy to plasma electrons and 
plasma ions, and the formalism we now develop will allow us to compute the electron-ion 
energy splitting in a systematic and unambiguous fashion. We take the plasma to have an 
electron temperature Tg = and a common temperature Tj — for all the ions, in 
which case the Fokker-Planck equation for the distribution / of an impurity species has the 
form 



d 



(2.1) 
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where v = p/w is the velocity of an impurity particle with momentum p, the explicit sum 
runs over all the particle species b in the background plasma, and the summation convention 
is used for repeated vector indices k and £. As we shall describe more fully, the diffusion 
coefficient C^^ has been analytically calculated to leading and next-to-leading orders in 
the plasma density in BPS [2] or more precisely, to orders In and g^ in the generic 
dimensionless plasma coupling constant g = c'^k/A'kT. We use rationalized electrostatic 
units, so that this parameter is the Coulomb energy of two particles of charge e a Debye 
distance 1/k apart divided by an average temperature T. 

With our conventions, the number of impurity particles is given by 

and their kinetic energy and momentum appear as 

and ,^ 

Since the right-hand side of the Fokkcr-Planck equation (2.1) contains an overall total mo- 
mentum derivative, it does not contribute to the time rate of change of the particle number 
— the Coulomb collisions in the plasma preserve particle number. When the electrons 
and ions arc at common temperature T = the terms in the final square brackets 

in the Fokkcr-Planck equation (2.1) annihilate a thermal Maxwell-Boltzmann distribution 
[/ oc cxp{— /3 p^/2m}] of impurity particles — a collection of particles in thermal equilib- 
rium is not altered by their collisions with a background plasma at the same temperature. 
However, for those cases in which the ions and electrons have different temperatures, the 
"injected impurity particles'' attain a non-thermal quasi-static distribution that will be de- 
scribed shortly. Eventually this quasi-static distribution will relax into a thermal distribution 
as the electron and ion components themselves thermally relax. As we shall see, however, 
the impurity distribution has interesting effects on temperature relaxation at intermediate 
times. 

The stopping power can be extracted from the Fokker-Planck equation by considering a 
single impurity particle at moving with the velocity Vp. The corresponding distribution 
function is given by /p(r, p, t) — (27r/i)^(5(r — rp)5(p — Pp), and one can easily check that this 
distribution indeed gives A?" = 1 as it should. Inserting this single particle distribution into 
Eq. (2.1) and performing a partial integration, it is easy to see that the rate of energy loss 
of the particle is given by 



dE _ y-^ 



v^C^'iPp). (2.5) 
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To make the sign of this expression clear, we emphasize that it gives the rate at which the 
particle loses energy to the plasma [it is the negative of the time derivative of Eq. (2.3)]. 
Hence the stopping power, which is the energy loss of the particle per unit distance traveled, 
appears as 

(2.6) 



dx 



Vp dt 



In a similar manner, we can find the rate of change of the momentum by substituting the 
single particle distribution into expression (2.4), thereby giving 



dt 



(2.7) 



As performed in BPS, by calculating dE/dt and dP^ /dt to leading and next-to- leading order, 
we can invert equations (2.5) and (2.7) to the same order to obtain the coefficients C^^ of 
the Fokker- Planck equation. 



B. Longitudinal and Transverse Components of the Diffusion Tensor 

As described in detail in BPS, the isotropy of the background thermal plasma allows one 
to decompose the diffusion tensor as 

(p) = A{v) ^ + B,{v) \ {S^^ - v'v') , (2.8) 

where v is the magnitude of the velocity, v — |v|, with the velocity direction given by 
ir — v/v. We often take the independent variable to be the energy E — \ mv"^ and, with a 
shght abuse of notation, we shall also write Ab — Ab{E) and = Bb{E). As a matter of 
completeness, the ^-coefficients are provided in Appendix A, and their various limits can 
be found in Appendix B. For a homogeneous and isotropic source of impurity particles, the 
case we shall consider, the B-coefficients do not enter, although their analytic forms can be 
found in BPS [2] if desired. 

Let us return to the stopping power (2.5) of a charged particle. Since the velocity tensor 
multiplying the B-contribution is transverse — its contraction with v'' or vanishes — the 
rate of energy loss (2.5) of a projectile becomes 

Ab, (2.9) 

where we have now omitted the p subscript. The respective energy losses to the ions and 
electrons are given by separating this formula into the ion contribution described by 

A^y^Ai, (2.10) 



dE 
'dt 



d 



V — 
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and the electron part governed by Ae, so that^ 



and 



with their sum giving 



~df 



dEe 
dt 



V — 



1 d 

1 d 



d^ 
'dt 



/3em dv^ 

dE, dE, 
+ 



(2.11) 



(2.12) 



(2.13) 



dt dt 

The rates dEi,/dt were rigorously computed in BPS to the leading g'^lng'^ and sub- leading 
orders, and these results were then used to determine Ab to these orders. In a similar 
fashion, BPS also computed the rate of momentum change dP^/dt of a projectile to these 
orders to determine the other coefficients Bb- In this way, a Fokker-Planck equation was 
determined to these orders in an unambiguous manner with no undetermined parameters.^ 
In particular, we should emphasize that our Fokker-Planck equation describes a particle's 
energy loss including orders g^ In g'^ and g^ with no ambiguity. 

Rather than tracking an individual charged particle slowing down in the plasma, it is much 
simpler — and equivalent — to examine an isotropic distribution of particles. When the 
impurity distribution is isotropic, / is a function the magnitude of the momentum p — |p| or 
equivalently, of the speed v or energy E. In such cases, a momentum derivative of / produces 
a factor of the velocity vector whose contraction with the velocity tensor multiplying the Bb 
coefficients vanishes. Hence in the isotropic case, the Fokker-Planck equation (2.1) reduces 
to 



d d ^ Ab 
. V > — 

rii rhr ■^—'6 777, 



1 -I- ^ 9 

^bmv dv 



f{E,t)^0. 



(2.14) 



To avoid notational clutter, we define the total ^-coefficient by 



A{E) ^ A{E) + Ae{E) , 



(2.15) 



^ As noted in BPS, to the order in g in which we are working, namely to leading (g^lng^) and next-to- 
leading (g^) order, only the kinetic energy of the stopping ion enters, and a meaningfully separation into 
electron and ion energy components can be made. This is because of the trivial fact that the kinetic 
energy is independent oi g — it is of order g^. In addition to this kinetic energy, the impurity particle has 
potential energy interactions with the ions in the background plasma. The change in these interaction 
energies associated with the motion of an impurity particle in a plasma cannot be separated into different 
parts that are associated with the ions and with the electrons. This is because this potential energy starts 
out at order g, and thus its evolution, which involves interactions akin to those involved in the kinetic 
energy dE/dx, is of order g^ (modulo possible logarithms), an order that is higher than that considered 
in this paper. Thus it should be emphasized that at higher orders in g, such clean separation into energies 
deposited into well-defined, separate ion and electrons components cannot be performed. 

^ See BPS [2] for a full discussion of the range of validity of the Fokker-Planck equation constructed in this 
fashion. 
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and the temperature-weighted ^-coefficient by 

{TA{E)) = Ti + Te Ae{E) . 

Thus, 



d__d_ 
dt dv 
Using the operator forms 

and 



AjE) ^ {TAjE)) ^ d 



f{E,t)^0. 



dv dv V dE 



d d d 
dv dv dE ' 



we may express Eq. (2.17) in the form 
d 1 d 



dt mv"^ dv 



v'A{E) + {TA{E)) 



V d 
m dv 



or 



dt mv dE 



A{E) + {TA{E)) A 



f{E,t)^0. 



(2.16) 

(2.17) 

(2.18) 
(2.19) 

(2.20) 
(2.21) 



C. Asymptotic Solution 

As we shall see, to use these results to obtain an unambiguous formulation of the frac- 
tions of the total energy deposited into the ions and electrons, we first need to compute the 
asymptotic distribution into which an initial swarm of test particles relaxes in the presence 
of a background plasma of differing electron and ion temperatures. This quasi-static distri- 
bution will be a function of E (or equivalently of p), which we express in terms of a function 
S{E) as 

U{E) ^ Af e-'^^\ (2.22) 
where we choose Af to normalize the distribution to unity. 



-S{E) 



(2.23) 



(27r?i)3 

The function S{E) is determined by inserting the structure (2.22) into Eq. (2.21) which 
gives 

w r ^ ^ 

(2.24) 



_d_ 
dE 



E !^A{E) + {TA{E)) e-'^^^ ^ . 



One solution of the second-order differential equation (2.24) is obtained by requiring that 
the quantity in curly braces operating on exp{— 5'(£')} vanishes: 



AiE) - {TA{E)) ^ = , 



(2.25) 
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and the solution can be obtained by a simple integration 

S(E;T,,T.) = l^E'^^^. (2.26) 

Here we have temporarily indicated the explicate dependence upon the electron and ion 
temperatures to emphasize that when the ions and electrons are at a common temperature 
T = Tj = Tg, this solution reduces to the Maxwell-Boltzmann distribution 

S{E;T,T) = |, (2.27) 

and consequently a swarm of test particles simply relaxes to the background plasma equi- 
librium distribution. For the equal temperature solution (2.27), a simple analytic Gaussian 
integration evaluates the normalization factor defined in Eq. (2.23) as 



. (2.28) 

Expression (2.26) is indeed the physical solution for S{E). This is because having the 
solution (2.26) in hand, it is a matter of simple quadratures to construct the second, linearly- 
independent solution for our second-order differential equation (2.24). It is not difficult to 
then confirm that this second solution is not normalizable, and so our first solution is the 
only physically relevant solution. We can also see that this is the desired solution since, for 
equal temperatures, it relaxes to a thermal Maxwellian distribution 

The Maxwell-Boltzmann distribution has an average energy of 3T/2. However, for the 
ions and electrons at different temperatures, the swarm of test particles relaxes to the average 
energy 

E = J^ fv^.^e^J-si^)]. (2.29) 
J (27rn)3 2m ^\ \2m J / ^ ^ 

In this case, numerical integrations are needed to evaluate the normalization constant Af 

and the average energy E. Figure 4 plots the average final energy E for an a particle in an 

equimolar DT plasma with an electron density Ue — 1.0 x lO^^cm"^. The figure displays E 

as a function of the ion temperature Ti for various electron temperatures Tg. 



III. FORMAL SOLUTION 

A. A Homogeneous and Isotropic Source 

We shall assume that the background plasma parameters, such as its density and tem- 
peratures, change very little over distances that are large in comparison with the stopping 
distance of the charged impurity particles, and that the plasma parameters also change very 
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FIG. 4: Average energy E to which the a particle relaxes as a function of the ion temperature 
Tj for various electron temperatures Tg. The ascending curves describing larger values of E have 
the increasing electron temperatures Tg = 10 , 30 , 50 , 100 keV. When = Tj = T then E = ^T. 
The background plasma is equimolar DT with electron number density rie = 1.0 x 10^^ cm^'^. 

little during the stopping time. Thus the plasma is treated as homogeneous and static. In 
addition, we assume that the sources of the impurity particles are distributed uniformly in 
space and that they emit the impurity particles isotropically with a definite energy Eq. For 
example, the fusion process in a homogeneous DT plasma produces a particles uniformly 
in space and isotropically in angle with an initial energy of Eq = 3.54 MeV. Thus, instead 
of considering the motion of a single test particle, we compute energy partitions and final 
states of charged particles emitted isotropically with a definite energy Eq from a uniform 
distribution of sources. This greatly simplifies the problem in that we can employ the homo- 
geneous Fokker-Planck Eq. (2.21) except that it is now modified to include a time-varying 
source of particles of energy Eq: 



The number and energy densities, n{t) and S{t), are simply given by removing the spa- 
tial volume integrations from the previous definitions (2.2) and (2.3). The inhomogeneous 
Fokker-Planck equation (3.1) gives the time variations of these quantities: 



di 



_2__d_ 

mv BE 



E A{E) + {TA{E)) — f{E, t) = 6 {E - Eq) s{t) . (3.1) 





6{E-Eo) s{t) 



s{t) 



\/2m>E( 



(3.2) 



(27rn) 
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and 

£{t) = Eoh{t)-J j^,vl^[ME)+Aem + [T.A{E)+TeAe{E)]^'^f{E,t). 

(3.3) 

When the impurity source s{t) is turned on and then attains a constant fixed value Sq, the 
number density n{t) eventually increases linearly in time, 

n{t) = //^'MO^^^^P LJ'''^''^ 

— rioo t + constant , (3.4) 

where 

noo = 2^A/^^. (3.5) 



B. Asymptotic Solution to the Inhomogeneous Problem 



We turn now to obtain the asymptotic solution to (3.1) satisfying the initial condition 
that there are no impurity particles in the distant past. 

As a first step in obtaining the asymptotic solution of the inhomogeneous Fokker-Planck 
equation (3.1), we set 

f(E, t) = exp {-S{E)/2} g{E, t) . (3.6) 

Multiplying the resulting Fokker-Planck equation by exp{5'(£')/2} on the left yields a sim- 
ilarity transformation that converts the (velocity ~ momentum) differential operator struc- 
ture in Eq. (3.1) into 



H = 



d_ . _ vA{E) 
dp'^ 2{TA(E)) 



{TA{E)) 



vA{E) 
9p ' 2{TA{E)) 



\ d 

V • 7^ + 



SO that the new Fokker-Planck equation now appears as 



^^^ + H'^g{E,t)^5{E-E,)e'^^'^y's{t). 



(3.7) 



(3.8) 



Incorporating the boundary condition that the solution vanishes initially, the inhomogeneous 
differential equation (3.8) has a formal solution: 

g{E, t) = f dt' e-^(*-*') 5 (E - Eo) e^(^°)/' sit') . (3.9) 

J — oo 

Because of the operator nature of the formal solution (3.9), it is convenient to view 
functions in momentum space as vectors in an abstract real vector space and define an inner 
product by 



(V',X) 



(27r/i) 



^^(p)x(p) ■ 



(3.10) 
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With obvious partial integrations, it is straightforward to verify that H considered as an 
operator on this function space is Hermitian with this definition of the inner product. 
In view of our previous work, it is easy to check that 

0(p) = A/-V2 exp {-S{E)/2} (3.11) 

now appears as a zero mode of the operator if, 

H(p = Q (3.12) 

that has unit normahzation, 

(0,0) = 1. (3.13) 

Except for this zero mode function, the remaining spectrum of H is positive. This is true 
because, for any function '^(p). 



d^p {TA{E)) 



\ d vA{E) 



^(p)^ >0, (3.14) 



dp 2{TA{E))_ 

since an examination of our results for the A coefficients shows that {TA{E)) > 0. The 
equality in Eq. (3.14) holds only if 



d ^ vA{E) 



V^(p) = 0. (3.15) 



dp 2{TA{E))_ 

The spherically symmetric solution ip{p) = i^{\p\) is clearly the previous zero mode function 
■0(p) = ((){E). Hence within the class of isotropic solutions — the only class that is relevant 
to our work — there are no other zero modes of H and all its other eigenvalues are positive. 
Since the operator H is Hermitian, 

0/7 = 0. (3.16) 
In view of this adjoint equation, it follows that 

0e-^(*-*') = l, (3.17) 

as one easily verifies by taking the time derivative. 

Except for this zero mode, we have shown that the other eigenvalues of the Hermitian 
operator H arc positive. This positivity constraint must be obeyed, for otherwise the Fokkcr- 
Planck would have diverging "runaway" solutions at large times. The operator that projects 
out the zero mode is obviously the outer product of the zero mode vector with itself, 

P = 0(g)0, (3.18) 

and we write the complement operator as 

Q^l-P, (3.19) 
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where the first term in (3.19) is the unit operator on the function space. By definition, the 
operator P acts on an arbitrary function as 

PV'(p) = (0 ® 0)v^(p) = 0(p) (0, ^) . (3.20) 

We now see that the unit operator in the form P -\-Q acting on g{E, t) in (3.9) produces 

g{E,t) = 0(p) j J^sMHE'-Eo) e^(W f_Jt's{t') 

+ f dt' e-^(*-*') Q5 (E - Eo) e^^^°^/^ s(t') . (3.21) 

J —00 

The momentum integral in the first term of (3.21) is easy to evaluate, 

As for the second term, since the operator Q selects out the positive eigenvalues of H, an 
integration by parts can be performed to produce 

r dt'e-"^'-''^Q5(E-Eo) e^^'''^^/^s{t') = ^ Q5 (E - Eo) e^^^°^/^ s(t) 
J-00 H 

- f dt'e-''^'-''^^QS{E-Eo)e^^^°^/^s{t'). 
J-00 ^ 

(3.23) 

We now assume that the source s{t) is adiabatically turned on and attains the constant 
value s{t) — So at late times. In the asymptotic limit, the rate s{t) is therefore vanishingly 
small and the second term in Eq. (3.23) may be neglected. We may also replace s{t) by its 
asymptotic value Sq in the first line. Hence, upon multiplying g{E,t) by exp{— 5'(£')/2} to 
return to the function f{E, t), we obtain 

fiE, t) = AAe-^(^) ^^gp dt' sit') + m , (3.24) 

where 

f{E) = e-^(^)/2 ^Q8{E- Eo) e^^^")/^ _ ^3 25) 

H 

Upon integrating (3.2), we can write this asymptotic late time solution more suggestively as 

f{E,t) = n{t)f^{E) + f{E), (3.26) 

with foo{E) — M e~^''^\ We emphasize that expression (3.26) is the asymptotic late-time so- 
lution to the inhomogeneous Fokker-Planck equation since the term involving the derivative 
s{t) is omitted. 
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FIG. 5: The waterfall analogy: the small blue rocks represent the plasma electrons while the larger red 
rocks are the plasma ions. The motion of the 'water' represents the evolution of the impurity ions that 
are injected into the background plasma. As 'water' falls down the electron-ion slope at a constant rate 
determined by f{E), energy is deposited into electrons and ions. At the bottom of the fall is a lake into 
which the excess 'water' drains and whose height n(t) rises linearly with time. 

C. Energy Deposition 

Before presenting an explicit version of the formal solution (3.26), we pause to describe 
its physical interpretation and its relation to the ways in which the stopping charged particle 
deposits its energy to the background plasma. At large times, the phase-space density has 
the time-independent contribution foo{E) = M exp{— S'(-E')} into which any set of initial 
test particles must relax, and the first term of (3.26) describes this distribution normalized 
to the correct density n{t). There remains a time-independent part f{E) that describes the 
stationary process of particles losing energy to the background electrons and ions as particles 
pass through "energy bins" from the initial energy Eq to the final asymptotic distribution. 
The situation described here can be pictured as the flow of water over a rocky waterfall that 
slows the motion of the water as it descends. The initial rate of flow of the river corresponds 
to the rate n{t)\ the height h of the waterfall giving a potential energy proportional to gh 
corresponds to the initial energy Eq. The energy dissipated in the fall corresponds to the 
energies lost to the ions and electrons. The final flow into a horizontal lake corresponds to 
the build up of the particles in their final distribution described by foo{E). This analogy is 
depicted in Fig. 5. 

1. Energy Splitting 

Upon inserting the form (3.26) into Eq. (3.3), we can identify the asymptotic constant 
rates of energy loss as 

En^ = [Eo-E,-E,]h^, (3.27) 
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in which 



and 



El 

Eq 



El 

En 



1 



1 



(277/1)3 



vAi{v) 



V Ae{v) 



d_ 
dE 



HE) , 



m . 



(3.28) 



(3.29) 



(27r/i)3 

Thus Eqs. (3.28) and (3.29) are the constant fractions of the original energy Eq deposited 
into ionic energy Ei and electronic energy Ee — the energy losses analogous to those of the 
water passing through the rocky waterfall. 



2. Plasma Heating and Energy Exchange 



exchange 

The original energy goes into energies lost to the ions and electrons, with the remainder 
the average final energy E of an impurity particle. For a background plasma with the 
ions and electrons at a common temperature T, E — 37/2, and the result (3.27) becomes 
obvious. 

When the electrons and ions have the same temperature T — — T^, the slowing 
down of fast particles in the plasma gives a steady-state heating rate per unit volume V — 
[El + E^ hoo- This heating raises the temperature T of the plasma, but in most cases, the 
rate of this heating is small in comparison with the slowing down time of the fast impurity 
particles, and so our quasi-steady-state computation is valid, with the temperature treated 
as a slowing varying function in our formulae. 

When the electrons and ions have different temperatures T^. and Ti, the situation may be 
quite different. In addition to the overall plasma heating V, the final ensemble of the impurity 
particles works to bring the electrons and ions to a common temperature. Returning to 
Eq. (3.3), we see that the final ensemble contribution produces energy density transfer rates 
to the ions and electrons given by 

d 



and 



Slit) 



+ 



vAj{v) 



vAe{v) 



{2'Khf 

Carrying out the energy derivatives yields 



1 + T, 



dE 



d_ 
'dE 



M e^v{-S{E)} n{t) , (3.30) 



M e-x^{-S{E)} n{t) . (3.31) 



(3.32) 



19 



and 



£e(t) = -{Te- T:) Q , (3.33) 

where 

C?. - = n{t) j ^ . ^^^^^^N exp {-S{E)} . (3.34) 

Since 

£,{t) + £,{t) = Q , (3.35) 

there is no net heating of the plasma. This process only brings the ions and electrons to a 
common temperature. 

In the absence of the impurity particle quasi-static equilibrium ensemble, the thermal 
relaxation rate coefficients are well approximated^ by 



Here 

(3.37) 

is the squared electron Debye wave number, and 

^1="^ (3.38) 

is the definition of the squared plasma frequency for particle a, with the electron squared 
plasma frequency ujI specified by a = e, while the total squared ionic plasma frequency uj'^ 
is the sum over all the ions in the plasma 

^? = E-^'- (3.39) 
In numerical terms, for an equimolar DT plasma, 

Cie = 3.13 X 10-2^ nl jln (^5.80 x 10^^ - l.Ssj cm-^ ps"^ , (3.40) 

in which the electron density Ue is measured in cm^^, the electron temperature in keV, 
and the overall units are (1.0 cm^"^)/(1.0 x 10~^^sec) as indicated. 

The total rate coefficient for electron-ion thermal relaxation is the sum Cje + Cj" . It is of 
interest to compare Cj" to Cig. Since Cj" is proportional to the number of impurity particles 



This is the sum of Eqs. (12.44) and (12.57) in BPS [2] as quoted in Eq. (12.12) except that a simple 
transcription error was made in the sum quoted in BPS in that the —7 — 2 in Eq. (12.12) should be 
replaced by —7 — 1. 
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FIG. 6: The ratio njCj^/nCie as a function of the electron temperature for an equimolar DT 
plasma with an electron density of 1.0 x lO^'* cm^^ for ion temperatures of 3 keV (blue), 10 keV 
(red), 30 keV (black), and 100 keV (green). 

that come into their final equilibrium state n{t), this comparison can be made independent of 
this density by evaluating the ratio of C^^/n{t) to Cie/rij. In Fig. 6 we plot this dimensionless 
ratio as a function of the electron temperature Te for various values of the ion temperature 
Ti ranging from 3 keV to 100 keV at an electron density rig = 1.0 x lO^^cm"^. Explicit 
calculation shows that the dependence of this ratio upon the electron density is weak. 
As rif. is increased from 1.0 x lO^'^cm"^ to 1.0 x lO^^cm"^, the greatest change in the ratio 
occurs for Ti ^ Tg : for Tj = 100 keV and Te = 3 keV, the ratio increases by 20%. 

We must add the caveat, already noted in the Introduction, that the discussion that we 
have just made applies only to the case in which the final alpha particle population is not 
large. Hence, although in some cases the ratios shown in Fig. 6 are of order one, the net 
effect of this new mechanism must be relatively small. 



D. Results in Terms of dE/dx 



The results (3.28), (3.29), (3.30), and (3.31) all have the generic structure 




(3.41) 



Here we may write 



(3.42) 
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use p = mv, and integrate by parts to obtain 



E - ! 



■fiE) 



m dv 



(3.43) 



As remarked in the Introduction [Eq. (1.6)], the expression to the right of f{E) in the 
integrand above is just v dEi^^^E) / dx . The integrand does not depend upon the direction of 
p. Thus the angular integration simply provides a factor of 47r. Using vdp = dE, we now 
have 

dE,,e{E) p2 



dE 



dx 2712^3 



f{E). 



(3.44) 



Thus all of our results involve a factor of the stopping power dE^^^/ dx for the ions or for the 
electrons, but the integration weight involves a more subtle function than those in the naive 
formulae (1.1) and (1.2) given in the Introduction. 



IV. EXPLICIT SOLUTION 



A. General Development 



We turn now to the explicit construction of the function f{E) from the formal expres- 
sion (3.25). We start by multiplying Eq. (3.25) by the (velocity ~ momentum) differential 
operator structure in Eq. (3.1). Passing this operator through the factor exp{— S'(-E)/2} is 
equivalent to the similarity transformation that converts it into the operator H. Hence, 



d_ 
dp 



■ v 



A+{TA{E))-.^ 
V op 



f{E) = e-^(^)/2 Q e+^(^)/2 5 {E - E^) so , (4.1) 



and remembering Eqs. (2.18), we see that this is equivalent to 



dv m 



A+{TA{E)) 



d_ 
'dE 



m 



(4.2) 



In the second equality we employed the definitions (3.18) and (3.19) of the operators P 
and Q and in the last line used the result (3.22). Obviously, a trivial first integral of this 
differential equation exists. Since the constant of integration must be chosen to make f{E) 
vanish at large E, this first integral reads 

(4.3) 



A{E) + {TA{E)) 



d_ 
'dE 
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where 9{x) is the unit step function that vanishes for x < 0. Note that, in view of the 
normahzation (2.23), 



and so the sum of the terms in the curly braces in Eq. (4.3) vanishes when E ^ 0. This 
is in accord with the fact that these terms on the right of Eq. (4.3) were produced by the 
integral of a derivative on the left-hand side of Eq. (4.2), a derivative of a quantity that 
vanishes at both E = and E = oo. Moreover, since the curly braces vanishes at E — 0, 
the right-hand side of Eq. (4.3) is finite at this end point as it must be. 

At this juncture, it is convenient to remember the definition (3.5) of rioo which can be 
expressed as 

mEo Ti'^h^ 



So = 



itiEq 



and to simplify the notation by writing 



so that we have 



Thus, Eq.(4.3) now reads: 



A{E) + {TA{E))-^ 



POO 

/ dE'VE^Jfe-^^^'^^l. 
Jo 



(4.5) 

(4.6) 
(4.7) 



m 



Eq 77,00 mEoE 



i^d{EQ-E) - J°° dE'^JTe-'^'"'^^ . 



To solve this differential equation, wc set 



-S{E) 



because then 



A{E) + {TA{E))^ 



d 



/(E)=e-^(^) {TA{E)) — -g{E)- 



(4.8) 



(4.9) 



(4.10) 



Since the integrating factor involves exp{-t-5'(i?)}, which exponentially increases without 
bound as the energy increases, to obtain a finite well-defined result we must integrate over 
the range E' — to E' — E and obtain 
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B. Energy Fractions and dE/dx 



The customary expressions for the energy fractions in terms of the stopping power 
[Eqs. (1.1) and (1.2)] emerge for low temperatures. To see this, we note that in this case, 
the energy integration range is very large on the scale of the temperature, and that the work 
of Appendix B shows that the electron contribution dominates over most of this range so 
that we may approximate 

5(E) (4.12) 

Moreover, the sum rule (4.7) imphes that the terms in the curly braces in Eq. (4.11) cancel 
when E' E, where £" is a lower energy limit that is on the order of the electron temperature 
Te. On the other hand, the integral in the curly braces in Eq. (4.11) is exponentially small 
when the integration variable E' is somewhat larger that the electron temperature T^. Hence, 
in the low temperature case, Eqs. (4.9) and (4.11) provide the approximate solution 

^^^^ - ^ f ^A'h - ^'^] v^<^' - - ^ 

(4.13) 



exp \-^{E- E')] = Te -4 exp (-1 {E - E') } , (4.14) 



Repeatedly using 

'""{~T, ^" "^j ^'dE' i Te 
and repeatedly integrating by parts, shows that the leading term in the low temperature 
case is given by the upper-limit contribution of the first term in this sequence: 

TT^h^ Te 1_ 

^ 'E {TA(E)) 

Placing this approximate result in the generic form (3.44) to evaluate the energy fractions 
(3.28) and (3.29) yields 

E, f^" dE Te dE, 



- ^ § (fJlEvA^^ - e)<E - ^) . (4.15) 



Eo Je Eo {TA{E)) dx 



{E) , (4.16) 



and 



E^ dE Te dE^(p^ 

In the low temperature case, the generic relation (1.6) gives dE/dx ~ A and so 



Te Te 



{TA{E)) Ti dE,/dx{E) + Te dEjdx{E) ' 
In the equal temperature case. 



(4.18) 



T 1 



{TA{E)) dE/dx{E) ' 

and, setting Eq — >■ 0, the low temperature expressions (4.16) and (4.17) reduce to the 
commonly used Eqs. (1.1) and (1.2) discussed in the Introduction. 
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C. Equal Electron and Ion Temperatures 



The case in which the ions and electrons have the same temperature, = T^. = T, is 
simple in several respects. First of all, it is physically simpler because the final distribution 
of the stopping charged particles is the Maxwell-Boltzmann thermal equilibrium distribution 
of the background plasma, 

exp{-S(i?)} = exp|-|| . (4.20) 

Thus, the energy transfer processes (3.30) and (3.31) do not appear because, with Eq. (4.20) 
holding, the combinations in the square brackets in these equations annihilate exp{— 5'(£')}. 
Thus, only the energy partitions Ei and Ee need to be examined, and these obey the obvious 
sum rule 

^T^Eo-E,-E,, (4.21) 

to which Eq. (3.27) reduces. Secondly, it is mathematically simpler because there is no need 
to find an explicit solution to Eq. (4.11) because Eq. (4.8) reduces to 



1 + tA 

dE 



Eorioo ° EA{E) mEo 

3/2 



Ktt^ ( ^) r dE' V2^e-^'^' . (4.22) 



EqEA{E) 

The operation in the square brackets that acts on f{E) on the left-hand side of this equation 
is just that which appears in the energy partitions (3.28) and (3.29). 

Placing this expression into the energy partitions Eqs. (3.28) and (3.29) and changing 
the momentum integration into an integration over energy expresses the fractional energy 
loss into ions and electrons as 

dE A,{E) f°° dE A^{E) 2 p^/^ 



and 



^ _ n dE Mm _ rdE ME)_ 2J^ r .r,,ry.-^E' u OA) 

Eo ~ Jo Eo A{E) i Eo A{E) ^ Je ' ^ ^ ^ 

where /3 = 1/T. Adding these equations gives 



[^i + ^e] = Eo- dE / dE'^e-^^' 

Jo VTT Je 

Eo-'^-^ rdE'E"/\-^^' 

Jo 



Eo-^T, (4.25) 
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where the second hne follows from a partial integration, and the last line from the definition of 
r(5/2). This is just the obvious result of energy conservation previously stated in Eq. (3.27). 

The results (4.23) and (4.24) can be simphfied for their exphcit evaluation. Writing these 
results with a trivial rearrangement of the terms presents them as: 



El 

Eq 



I 



^° dE AjE) 
Eo A{E) 



1 - 



2^3/2 



TT 



dE'^'e-^^' 



L Eo A{E) Je 



(4.26) 



and 



E. 



Eo 



dE AejE) 
Eo Jo Eo A{E) 



1 - 



2^3/2 



TT 



dE'VWe-^^' 



dE AejE) 2/3^/2 
E^ A{E) Je 



dE'VE^e-^^' . 



(4.27) 



As we shall see, the second set of double integrals in Eqs. (4.26) and (4.27) are exponentially 
small. Hence it suffices to use the simple bounds 



ME) 



and similarly 



AjE) 

A{E) A^{E)+Ae{E) 

Ae{E) 



<1, 



A{E) 



< 1. 



(4.28) 



(4.29) 



Using these bounds, we encounter 

dE 2/33/2 



Eo Eo y/n Je 



dE'VE^e-^^' = 



2^3/2 



2^3/2 



TT 



Eo 
00 

Eo 



dE'VE^e-^^' I ^ 

Eo ^0 

E' 



E' 



dE 



dE'y/E'e-^^' 



E, 



1 







. (4.30) 



The variable change E' — Eo{x + l) presents this as 



£;3/2 r dx{l + xy/'xe- 

V'^ Jo 



■0Eox 



-pEo 



(4.31) 



with the evaluation on the right-hand side following from the fact that (5Eo ^ 1 so that only 
small X regions contribute justifying the replacement (1 + x)^/^ — > 1. Hence we indeed find 
that the additional double integrals in the energy fractions (4.26) and (4.27) are exponentially 
small, and so with very good accuracy we may write these fractions as 



El 

Eo 



Eo 



dE A.jE) 
'Wo A{E) 



1 - 



2^ 



dE'VWe-^''' 



E 



(4.32) 
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and 



En 



I 



^° dE AejE) 
Eo A(E) 



/ dE'^/We-''^' 



Since 



we may write 



2^2 



6/2 poo 

=^ / dE'VWe-"''' ^1, 

71" Jo 



1 - 



2^3/2 



r dE'y/E' e-^^' = ^ ^ dE' e"^^' , 
Je Vtt io 



^TT Je 

while partial integration gives 

2 ^3/2 ^i? 



A io 

Hence, using the definition 



dE' ^/E'e-^''' ^-xl^ e-^^ + J ^ P dE' ^ e'^^' 



TT 



'E' 



2 2 

erf(a;) = / dy 
VTT Jo 



of the error function, we may write the results (4.32) and (4.33) as 



Eo 



dE A,{E) 



Eo Jo Eo A{E) 



crf(V^£;)-A/^e-'3^ 



TT 



and 



^0 dE Ae{E) 



erf(V^£;)-A/^e-^^ 



TT 



£;o Jo ^0 A{E) 
These are the results quoted in Eqs. (1.7) and (1.8) in the Introduction. 



(4.33) 



(4.34) 



(4.35) 



(4.36) 



(4.37) 



(4.38) 



(4.39) 



D. Differing Electron and Ion Temperatures 

As we have seen, when the ion and electron temperatures of the background plasma 
differ, Tj ^ Tg, both the physical interpretation is richer and the mathematics becomes more 
difficult. With different temperatures, there is the additional physical process in which the 
final distribution of stopped injected impurity particles works to bring the electrons and ions 
into thermal equilibrium at a common temperature T — Ti = Tg. Moreover, mathematically, 
we must now work with Eq. (4.11). 

We use Eq. (4.11) to return to the f{E) function, and insert the result for f{E) into 
Eqs. (3.28) and (3.29) to compute E^/Eq and Eg/Eo- To simphfy the resulting formulae. 
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and place them in a form that parallels those for the previous equal ion-electron temperature 
case we note that 



AiE) 

Hence, with the definition 

G(Ti, Te; Eq 



dE 



-S{E) 



+ 



{TA(E)) 



-S{E") 



(4.41) 



the energy loss fractions may be expressed as 



E^ 

Eq 



Te-T, 

En 



G{Tj,Tf.; Eq) 



+ 



(4.42) 



and 



Ee 
Eq 



+ 



Eo {TA{E 



GiT„T,;Eo) 



Eq 

°° dE Te Ae{E) 



(4.43) 



The second fines in the results (4.42) and (4.43) are straightforward generalizations of 
the common ion and electron temperature forms (4.23) and (4.24). The first lines of the 
new results (4.42) and (4.43) cancel when they are summed, so that 

pEq poo poo 

E, + Ee ^ + dE-Jf dE dE' e-'^^'^'K (4.44) 

Jo Jo J E 

Upon interchanging integrals, 

poo roo roo pE 

dE dE'VWe-^^^'^ = / dE^Ee-^^^^ / dE' 

Jo J E Jo Jo 



f 

Jo 



dEEVEe-^^^K 



(4.45) 



Hence, on passing from an integration over energy to an equivalent momentum integral and 
reverting to the corresponding normalization factor J\f, we have 



El + Eg — Eq — 



(27r/i)3 



(4.46) 
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or, in view of Eq. (2.29), 



El + Ef, — Eq — E , 



(4.47) 



in which E is the average energy to which an impurity particle relaxes. This result is in 
accord with the previous Eq. (3.27). 

The final energy integrals in Eqs. (4.42) and (4.43) run from E — to E ^ oo. In each 
case, the final integration region involves the exponentially small factor exp{—S{Eo)} ~ 
exp{—Eo/T}, where T is a typical plasma temperature. This is a very small factor, and 
hence this upper portion of the integration region may be safely neglected to write the results 
as 



El 

Eq 




GiT„n;Eo)+ ^ 
Jo ^0 



dE TiAi{E) 



1 - 



G(Ti, Te; Eq) 



{TA{E')) 

^° dE T\A{E) '■^ 
^ {TA{E 



% [ dE'^'JIe-^^'''\ 
)) Jo 



(4.48) 



and 



Ee 
Eq 




I, Te] Eq) + / 

Jo 



Q Eq {TA{E')) 

^°dETeAe{E) 



Eq {TA{E')) Jq 



dE'VE^Jfe-^^'''^ 



(4.49) 



Here we have invoked the sum rule (4.7) to write the second equalities above. 

The work in Appendix C shows that the function G can be approximated, with an 
accuracy of a few percent, by 



GiZ,n,EQ)-J^ dEE^^^^j^e -WWAmJo ^' 

, A(^o)A(^o) 



-S{E") 



A^{Eq 



(4.50) 



Since the integration in the first line is over the finite interval (0, £"0) and since it involves 
only nested integrals, rather than a three-dimensional integral with an arbitrary integrand 
that involves an general function of three variables, its numerical evaluation is not difficult. 

The explicit forms for the A coefficients reviewed in appendix A now enable the explicit 
computation of the energy ratios E^/ Eq that are presented in Fig. 3 and the tables of 
Appendix V. 
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V. SUMMARY AND CONCLUSION 



Wc have developed a formalism that enables the calculations of the energy fractions that 
a fast particle deposits to the ions and electrons when it slows down in a plasma of ions and 
electrons that have different temperatures. Such calculations have not be done previously. 
Our work applies to background plasmas that are weakly to moderately coupled — the range 
of validity of this restriction was discussed in the Introduction. 

Since the background plasma is not in thermal equilibrium, a fast particle ends in a 
"schizophrenic" distribution which we explicit compute in Sec. IIIB. As described in 
Sec. Ill C 2, the final non-thermal distribution of the initial fast particles provides a mech- 
anism to bring the differing electron and ion temperatures to a final common temperature, 
a process that now appears in addition to the usual electron-ion relaxation interaction. 

Although our general method applies to the slowing of any fast particle in an arbitrary 
background plasma, we are specifically interested in DT nuclear fusion, and thus we present 
explicit numerical results for an initial 3.54 Mev alpha slowing in an equimolar DT plasma. 

For the case of equal ion and electron temperatures, the energy fractions Ei/Eq and Ee/Eo 
that wc compute arc in agreement with previous work to leading accuracy, but our results 
are more precise because we also compute the exact coefficient of the first non-leading term 
which is proportional to the plasma density n. The comparison between our and previous 
results for the equal temperature case was discussed in the Introduction and shown there in 
Figs. 1 and 2. 

In order to motivate and give the fiavor of our results for the general case in which the ions 
and electrons in the background plasma have different temperatures. Fig. 3 was presented in 
the Introduction. The table that follows gives detailed results for the energy fractions E^/Eq 
and E^/Eq for an alpha particle with an initial energy of 3.54 MeV slowing in equimolar DT 
plasmas of three different densities and a variety of temperatures. 
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DT FUSION ALPHA PARTICLE ENERGY DEPOSITED 
INTO THE IONS FOR VARIOUS PLASMA CONDITIONS 



TABLE I: Ei/Eq for a density n,, = 10^^ cm over a range 

of electron and ion temperatures that are measured in kcV. 





J-e = 10 


on 

Je = 20 


on 

Je = 30 


Je = 50 


Je = J^OO 


ie = iOO 


10 


0.248 


0.404 


0.513 


0.660 


0.834 


0.936 


30 


0.234 


0.389 


0.497 


0.644 


0.821 


0.930 


50 


0.220 


0.374 


0.481 


0.628 


0.807 


0.922 


100 


0.185 


0.336 


0.440 


0.586 


0.769 


0.892 


200 


0.126 


0.263 


0.361 


0.502 


0.689 


0.827 


300 


0.079 


0.197 


0.285 


0.418 


0.607 


0.760 



TABLE II: Ei/Eq for a density Ue = 10^^ cm ^ over a range 

of electron and ion temperatures that are measured in keV. 





re = io 


Te = 20 


Te = 30 


Te = 50 


Te = 100 


Te = 200 


10 


0.267 


0.421 


0.531 


0.675 


0.843 


0.939 


30 


0.252 


0.406 


0.515 


0.659 


0.830 


0.933 


50 


0.236 


0.391 


0.499 


0.643 


0.816 


0.925 


100 


0.200 


0.352 


0.457 


0.601 


0.778 


0.895 


200 


0.139 


0.278 


0.376 


0.516 


0.698 


0.831 


300 


0.089 


0.209 


0.299 


0.431 


0.617 


0.764 



TABLE III: Ei/Eq for a density He = 10^^ cm ^ over a range 

of electron and ion temperatures that are measured in keV. 



Ti 


Te = 10 


Te = 20 


Te = 30 


Te = 50 


Te = 100 


Te = 200 


10 


0.293 


0.446 


0.555 


0.694 


0.854 


0.942 


30 


0.276 


0.430 


0.539 


0.679 


0.841 


0.936 


50 


0.260 


0.415 


0.523 


0.663 


0.827 


0.928 


100 


0.223 


0.375 


0.481 


0.620 


0.789 


0.899 


200 


0.159 


0.299 


0.398 


0.534 


0.710 


0.836 


300 


0.106 


0.228 


0.318 


0.449 


0.630 


0.770 
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Appendix A: The A-CoefRcients 



The Fokker-Planck equation described in the text involves two scalar coefficient functions 
with only one of them, the A coefficient, entering into our problem of the partition of 
the energy loss of a fast charged particle into the ions and electrons in the plasma. The 
Fokker-Planck equation, and the coefficients Ai and Ae coming from the ions and electrons 
that are needed for our problem, were discussed extensively in BPS [2]. There a method 
of dimensional continuation was employed to compute the Ab which enables the short- 
distance, point Coulomb scattering to be joined with the long-distance, collective force in an 
unambiguous fashion that has no double counting. This method was used to evaluate the 
Ab both to leading and to subleading order — roughly speaking — to order n In n as well as 
n, where n is the plasma number density (made dimensionless by the adduction of suitable 
parameters). For completeness, we present here the results of BPS. Since their derivation is 
subtle, it cannot be sketched here. 

The coefficient for the interaction of an "impurity particle" of energy E or velocity Vp, 
{E = mpVp/2) with the species h of the background plasma may conveniently be written as 

Ab{vp)=At{vp)+At'^{vp), (Al) 

which is the same as Eq. (10.25) of BPS, with 

^>p) = ^bM) + A"r(^p) , (A2) 

which is the same as Eq. (9.6) of BPS. Here Ab{vp) has two terms. The first accounts for 
the hard Coulomb scattering in the classical hmit, while the second accounts for the collec- 
tive, long-distance effects, which are entirely classical. The term A^'^{vp) is the quantum- 
mechanical correction to the scattering that vanishes in the limit in which Planck's constant 
vanishes, — > 0. 

The first classical piece is given by 

A,sM = ) exp|--/3(,mfeVpM 



-ln(/?.^ir^-^)-27 + 2 
47r nipb 1 — u 



(A3) 



which is contained in Eq. (9.5) of BPS. The reduced mass nipb of the projectile (p) and 
plasma particle (b) is defined by 

' ' (A4) 



nipb rup nib 
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The second part of the classical contribution is given by 

=2 • rl 



"'^^ 47r27ry_i PtotaiKcos^) ' \ i 

which is contained in Eq. (7.26) of BPS. Here Ptotsiiy) is the spectral weight, 

Ptotal(^^) = ^,Ph{v) , (A6) 



with 



Pb{v) = 4 exp|-^/3bm6V^| , (A7) 

as introduced in BPS Eqs. (7.9) and (7.10). With these definitions, it is not hard to show 
(as is explicitly done in BPS) that the sum + A^^ is independent of the arbitrary wave 
number K that was introduced for computational convenience. The function F{vpC0s9) is 
related to the classical dielectric function e{k, kvpCOsO) by 

e{k, kvp cos 6) ^ k"^ + F{vp cos 6) . (A8) 

Here, consistent with our leading orders evaluation, the dielectric function corresponds to 
the classical limit of the quantum ring sum. Hence the complex- valued function F{v) is 
defined by 

F{v)^- Tdu-P^^^^^. (A9) 

Equations (A8) and (A9) are the formulae (7.7) and (7.8) of BPS. 

The quantum correction is contained in Eq. (10.27) of BPS, and it reads 

exp{4Am.(.,-.,.f}(l-^-i— ) 



PbmbVpVpb 



+ exp <! -]r/3bmb {vp + Vpbf 



2 J V PbmbVpVpb 



(AlO) 



Here ip{z) — dlnr{z)/dz and, with the rationalized Gaussian units that were used by BPS 
(and which we continue to use) where the Coulomb potential energy between charges Ca and 
Cb a distance rab apart is given by y = eaeb/{4:7rrab), the formula contains the dimensionless 
quantum coupling 
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The following figures illustrate the behavior of the ^-coefficients for an equimolar DT 
plasma with an alpha particle projectile of kinetic energy E. Figures 7 and 8 plot the electron 
and ion components Ae, Ai and their sum A = Ae + Ai for a plasma with electron number 
density Ug = 1.0 x lO^^cm"^, electron temperature Tg = lOkeV, and ion temperatures of 
Ti = lOkeV and 100 keV, respectively. 

0.15 I 




E (MeV) 



FIG. 7: The coefficients Ai{E) (dashed blue), Ae{E) (dotted green) and A{E) (sohd red) as functions of 
the kinetic energy of an a particle projectile. The background plasma is equimolar DT with electron 
density rie = 1.0 x 10^^ cm~'^ and electron-ion temperatures = lOkeV and = lOkeV. 



0.05 




12 3 

E (MeV) 

FIG. 8: As in Fig. 7, except Tg = lOkeV and = 100 keV. The crossover energy E ^ Ec where 
Ae{E) = Ai{E) is about the same in both Figures; however, the peak value of the coefficient Ai is inversely 
proportional to Tj. 

Figures 9 and 10 illustrate the number density scaling of the ^-coefficients by plotting 
Ae{E)/ne and Ai{E)/ne, as a function of the a particle energy E, over a wide range of 
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electron densities: rig = 10^^, 10^^, and 10^^ cm ^. As before, the electron temperature is 
Te = lOkeV and the ion temperatures are Tj = lOkeV and Tj = 100 keV, respectively. 
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E (MeV) 

FIG. 9: The ^-coefficients for electrons and ions as a function of the a particle projectile energy E in 
an equimolar DT plasma with equal electron and ion temperatures, Tg = Ti = 10 keV. The solid lines 
correspond to Ue = 1.0 x lO^^cm"'^, the dashed lines to rie — 1-0 x 10^^ cm"'^, and the dotted lines to 
fie = 1.0 X lO^^cm"'^. In each case, the ^-coefficient has been rescaled by the corresponding number density 
Ue- The slowly rising (red) curves are those for Ae, while the sharply peaked (blue) curves are for Ai- 




E (MeV) 



FIG. 10: As in Fig. 9, except with Te = lOkeV and T, = 100 keV. 



Because the v4-coefficients are proportional to the Debye wave number squared, a quantity 
proportional to He, it is no surprise that Ai and Ae approximately scale with Ue- The Debye 
wave number also appears inside the logarithm and the dielectric function, and for electrons 



35 



this produces a much more pronounced effect than for the much heavier ions: while Ai/ue 
is almost independent of rig, the electron component Ae/fie varies by a factor of two over 
the range of Ue- 



Appendix B: Asymptotic Limits 

We shall extract the large and small energy limits of the Ah{vp) function for the various 
plasma species b from the general expressions in BPS [2]. The energy is given by £■ = 
where rup and Vp are the mass and speed of the particle moving through the plasma, the 
projectile p. We shall obtain the large and small limits of the projectile energy E as compared 
to a typical plasma temperature T. 



1. E ^ T: Electrons and Ions 



In the low velocity limit, Ab{vp) vanishes linearly with Vp, and so we write 



■Up ^ : 



A.,.p,^^^(^\'\p[A,,At^], 



(Bl) 



with two constants and Af''^ . These two constants arise from the low velocity limit of the 
classical and quantum pieces of Eq. (Al). The classical piece has already been calculated by 
BPS, where it is contained in their Eq. (9.9), so there is no need to do it here. The result is: 



Ai 



In 



IGtt rupb 



epCbPbKo rub 



1-2, 

2 ' 



(B2) 



in which rUpb is the reduced mass defined in Eq. (A4) in the previous Appendix, and 7 = 
0.577 ... is Euler's constant, and 



«D ^^^>^l^ Xlfe ^b ■ 

To bring out the size of this classical part, we define a plasma coupling by 

epCb/SbK' 



9pb — 



An 



47rAn Tb ' 



(B3) 



(B4) 



in which Ad = 1/kd is the Debye length and Tb — 1/ f3b is the temperature of plasma species 
b. Then we may write 



Al 



ln(±^)-i-27 
Qpb nib J 2 



(B5) 
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which shows that > 0, since Qpb must be small for our perturbative computation to hold. 
Note that when the electron and ion temperatures are not vastly different, the ions dominate 
in the low velocity limit (Bl) by a factor (mj/me)^'^^. Moreover, since ~ 1/T, this ionic 
contribution to the ^-coefficient has the temperature factor and thus increases as the 

ion temperature is lowered. The corrections to the low energy limit (Bl) are of relative order 
E/T. 

The low velocity limit of the quantum correction Eq. (AlO) above was not previously 
calculated in BPS because there the low velocity limit of dE/dx was used only to compare 
with a computer simulation involving classical dynamics, and therefore the quantum correc- 
tion was not needed. The needed quantum part is contained in Eq. (10.27) of BPS which 
provides the limit 



^ 



— /3fe"^6 ^ dvpb Vpb exp I - ^ l3bmbvli, | [2Reip {1 + irjpb) - In 77^^,] . 

(B6) 



To bring out the character of Eq. (B6), we introduce a thermal velocity Vb by 

1 -2 3^ 
2'^bVb = 2 



or 



PbTUb ' 
CpCb 

Airhvb 

We then change the integration variable, 

CpCb ^p^b 



and a corresponding quantum parameter 



Vpb 



Vpb 



u^ripbVbU. 



to obtain 



^ iVpb) = -vlb I duuexp\~r]l^u^ 



(B7) 
(B8) 

(B9) 
(BIO) 
(Bll) 



If we introduce the Bohr radius Oq = Airh'^/e^'me and use the average squared thermal 
velocity definition (B7), we can write 



_2 1 fepeb\'^ rrib 1 1 /'CpCbX'^ rrib 27eV 

^'^^ 3 V / rrie Tb Atiuq 3 V / rrie Tb 



(B12) 
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Thus for the charge and mass of a typical projectile particle such as an alpha particle and 
for a typical hot plasma, we see that for the electrons in the plasma fjp^ <^ 1, while for the 
ions in the plasma jy^j ^ 1 unless the ion temperature is somewhat larger than 10 keV. 

For fjp^ <^ 1, the exponential does not rapidly damp large u values, and so the relevant 
piece of the integrand is that with 1 where 



( 1 + - ) -^'(l) = -7, 



(B13) 



leading to 



_2 

~ -rj 



pe 



duu exp 



^ -2 2 



[-2 7 + lnw2] 



(B14) 



Adding this result to the classical limit (B2) gives the complete plasma electron contribution 
for a low energy projectile: 



^ « T T]l « 1 : 



Ae{Vp) = 



4tt V 27r 



1/2 



7 



(B15) 



For Tjpi ^ 1, the exponential rapidly damps large u values, and so the relevant piece of 
the integrand is that with u <^ 1 where 



2Reij 1 1 + - J +lnM' 



6 



(B16) 



and thus 



ti » 1 



2 2 



1 

27 



(BIT) 



Since this is a very small correction to ^4? > 0, it may be neglected, and we may use the 
pure classical limit (B2) for the ion contribution: 



Ai{vp) 



^p ^1. 



I3im, 



1/2 



2v. 



In 



167r m. 



pi 



27 

2 ' 



47r V 27r y 3 
The total contribution of the ions in the plasma in this case is obviously 



(B18) 



£;«r, f^J,»i: 



(B19) 



38 



0.02 0.04 0.06 0.08 0.1 

Ve (VMeV) 

FIG. 11: The ion contribution (green) plotted with the corresponding low-energy approximate form (red- 
linear) given by Eqs. (B18) and (B19). The plasma is equimolar DT with Tg Tj = lOkeV and = 
1.0 X 10^^ cm""^, and the projectile is an a particle. For these parameters, the plasma coupling is ge = 0.0006. 
Since the leading-order small-energy behavior is proportional to Vp, the graph of Ai against ^/E is linear in 
this region. 




0.5 1 1.5 Vs. 5 

Ve (VMeV) 



FIG. 12: The electron contribution (blue) plotted with the corresponding low-energy approximate form 
(red-linear) given by Eq. (B15). The plasma is the same as in the previous figure. Note that the linear 
approximation holds well into the DT fusion production energy of 3.54 MeV for the a particles. 

In Figs. 11 and 12 we plot the ion and electron ^-coefficients for an equimolar DT plasma 
with an electron density ne = 1.0 x 10^^ cm~^ and equal electron and ion temperatures 
Te = Tj = lOkeV against the square root of the projectile energy ^/E. We make this 
choice because in the small-energy regime the coefficients are linear in the projectile velocity; 
therefore, the graphs exhibit linear behavior until they start to depart from the low energy 
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limit. Figures 13 and 14 plot the ^-coefficients for an equimolar DT plasma with Te = 10 keV 
and Ti = 100 keV and with an electron density Ue = 1-0 x lO^^cm"^. 




0.05 0.1 0.15 0.2 0.25 

Ve (VMeV) 



FIG. 13: The ion contribution (green) plotted with the corresponding low-energy approximate form (red- 
linear) given by Eqs. (B18) and (B19). The plasma is equimolar DT with Tg = lOkeV, Ti = 100 keV and 
He = 1.0 X 10^^ cm~^, and the projectile is an a particle. 



> 
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FIG. 14: The electron contribution (blue) plotted with the corresponding low-energy approximate (red- 
linear) form (B15). The plasma is the same as in the previous figure. 

In Figs. 15 and 16 the temperatures are changed to Tg = 100 keV and = lOkeV. 
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0.15 ^ 




0.02 0.04 0.06 

Ve (VMeV) 



FIG. 15: The ion contribution (green) plotted with the corresponding low-energy approximate (red-linear) 
form (B18) and (B19). The plasma is equimolar DT with Te = 100 keV, T, = lOkeV and = 1.0 x 
lO'^^ cm^"^, and the projectile is an a particle. 



0.0004 
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FIG. 16: The electron contribution (blue) plotted with the corresponding low-energy approximate (red- 
linear) form (B15). The plasma is equimolar DT with Tg = 100 keV, T, — lOkeV and rie — 1.0 x 10^^ cm^^, 
and the projectile is an a particle. 

In Fig. 17 the coefficients Ae and Ai are plotted together for three different temperatures. 
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0.5 1 1.5 

Ve (VMeV) 

FIG. 17: The coefficients Ae and Ai are plotted together for three different temperatures for an equimolar 
DT plasma with an electron number density rig — 1.0 x lO^^cm"^. The ion contributions Ai peak to the 
left in the figure. The temperatures are: (i) Tg = lOkeV and T, = lOkeV (blue), (ii) = lOkeV and 
Ti = 100 keV (red), (iii) = 100 keV and T, = lOkeV (black). The electron contributions A^ for cases (i) 
and (ii) are almost equal, whereas for case (iii) Ae is very small. 



2. E^T: Total Ionic Contribution 



For the total ionic contribution, it is convenient to first work out the regular part of 
the long-distance, dielectric contribution because it is the same for both cases of classical 
and quantum-mechanical scattering. With a trivial integration variable change, Eq. (A5) 
presents this contribution as 

A<{v) = i^^ r dv V ^'^"^ F(v) In (^] (B20) 



where we now write 



so that 



a(^) = $^.a(^), (B21) 



PtotA^) = Pe{v) ^ pIv) , (B22) 

with the weight functions p given by Eq. (A7). Assuming that the charges of the ions do 
not differ greatly from the charge of the electron, then k1/ ti^ ~ Tj/Tg and the integrand of 
Eq. (B20) involves a factor that has the behavior 

^'(^) - ^ \ , (B23) 

Ptotai(w) 1 + pe{y)/pi{y) 1 + {meTj^/m.T^y/'^ exp {m,v'^/2T,} ' 

where rrii is a typical ion mass. Thus, defining a typical ionic thermal velocity v^^ by 

mi vl = T,, (B24) 
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this factor remains unity up to the critical velocity Vcrit defined by 



-crit = -^ln(^') , (B25) 



meT3 

after which it falls fairly rapidly to zero. The logarithmic factor in Eq. (B25) is typically 
about a factor of 10. So Vcnt is somewhat larger than an ion thermal velocity yet it is 
considerably smaller than the electron thermal velocity. 

In this region in which the factor Pi('y)/ptotai(^) of the integrand is non- vanishing, the 
function [Eq. (A9) above] 



Fiv) = - / du ^^"^^'^-^ (B26) 
has the form 

F{v) = F{v) ^kI + F,{v) , (B27) 

where 

F,{v) = - du ^'^""^ . (B28) 

This is so because the velocity which must be less than Vcriti is much less than the electron 
thermal velocity. Hence the electron part of F{v) takes on its low velocity limit, the electron 
Debye wave number squared k^. We place the form (B27) into Eq. (B20) to obtain 

Here we have replaced the arbitrary intermediate wave number K by the electron Debye 
wave number Ke because now 

^; ^ oo : ^ 1 _ ^ , (B30) 

Kg Kg 

where ^ 



and so ln(F(^;)/Kg) vanishes for large v. 
In order of magnitude. 



2 2 - 2 ^ 2 • (B32) 



Hence, since v?^ is much less than fcriti unless Tg is considerably larger than Tj, the final 
factor in the integral (B29), ln(F(f vanishes before Pi{v) / ptotaxiy) departs significantly 

from unity. Hence we simply take pi{y) / piota\{y) - 



el 1 



V 




^tM) ^^-2^1 In I I . (B33) 
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The discussion above shows that when Vp > t'crit, the hmits of the integration may be 
replaced by ±00. RecaUing the definition (B25) of the critical velocity Vcrit, and assuming 
that the projectile mass nip is about the same as the typical ion mass mj in the plasma, we 
can now state that 

Wc should note the convergence of the integral requires that the integration limits are to be 
taken in a rigorously symmetrical fashion with the integral performed between exactly —Vp 
and +Vp and then f p — )> 00 taken. It is now a simple matter to evaluate this limiting form. 
Adding a semicircle in the upper half plane of radius Vp gives a closed contour integral with 
no interior singularities that accordingly vanishes. Hence the value of the original integral is 
the negative of the integral over this large semicircle, an integral that is trivially performed 
using the limiting forms listed before. Thus 

^ > ^ 1^ (^] ■ = -^^,- (B35) 

With the long-distance, dielectric ionic contribution evaluated in the projectile high- 
energy limit, we can now compute the complete function Aiivp) in this limit. To do so, we 
must distinguish two cases for the remaining hard scattering contribution. 



^ > T In 



a. E^T, ri^- » 1 

As shown in detail in Sec. 10 of BPS, the classical scattering contribution dominates 
when the Coulomb parameter rjpi is large, with the first quantum-mechanical correction of 
relative order 

o _ (AnhVpV f V 2E 



pi 
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(B36) 



where a ~ 1/137 is the fine structure constant. In this classical limit, the scattering con- 
tribution is given by Eq. (A3). For the previous evaluation of the dielectric contribution to 
hold, we must choose K — so that this formula reads 
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Since miv'^/2 E ^ Ti, only small u values are significant. Hence we can approximate 
1 — u = 1 within the logarithm and extend the integration limit to w — >■ oo. With the 
variable change {^imiv'i/2) u—s^,we obtain the high-energy limit 



el CO? 4 



ds exp {— s^} 
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where we have used K^/Pirrii = ou"^. Here we have the integrals 

4 

and 
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Whence, 
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which, summed over all the ions in the plasma and combined with the previous long-distance 
result (B35) yields the total contribution from the ions in the plasma: 

£;»r, r?2. » 1 : 
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See Fig. 18 for a comparison of Ai with its asymptotic form at large energy. 



In this case, we have the limit 



Vp ^ 00 : 



which is contained in Eq. (10.42) of BPS. Adding this result to Eq. (B35) now provides the 
complete Vp ^ 00 limit for the ion part of the Ai coefficient: 
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FIG. 18: Ai (blue curve) vs its large-energy (red curve) asymptotic form. The plasma is equimolar DT 
with He = 1.0 X lO^^cm-^ and T^=T, = lOkeV. 

c. Rough Estimate For Either Case 

The overall factor in the contribution of the ions to either of the high velocity limits (B42) 
or (B44) of the ^-coefficient has the same form as that of the electron contribution (B59) 
given below, but with the major difference that the squared electron plasma frequency ujI 
is replaced by a sum of squared ion plasma frequencies uf, which are much smaller than 
the electron contribution by the ratio rrie/mi. For a rough estimate of the size of the ionic 
contribution for high energy projectiles, we approximate the logarithm in either limit (B42) 
or (B44) by a constant L of order one, and approximate uj^ ~ T^Ki/rrii and rrii ~ nip to 
obtain 



3. T <^ E <^ iripT/me: Electronic Contribution 

There is an intermediate range of projectile energies in which the projectile energy is much 
larger that the temperature, E ^ but yet not so large that we have E <^ {nip/ me) T ~ 
10^ T. We examine this range here. 

We again need to work out its long-distance, dielectric contribution, and its short-distance 
scattering contribution. 




(B45) 
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a. Dielectric Part 



In the energy range specified, the typical velocity in the dielectric function is small in 
comparison with the electron average thermal velocity and large in comparison with an ion 
average thermal velocity. Hence, in this range 

F{v) c^nl--L+m pt^Uv) . (B46) 



Here, in the dominant integration range. 



I ? T 

' - «1, (B47) 



and so we may simply write 

~«;2 + 7riptotai(^^). (B48) 

Moreover, in the dominant integration range, the imaginary part tt Ptotai(f ) is small in com- 
parison to Kg. Writing Eq. (A5) as 

2 • /•! 
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A^Jvp) ~ — / dcos^cos^ 



F{vp COS — F{—Vp COS ( 
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+ 



and using Eq. (B48) with the imaginary part treated to first order, 

=2 /•! r / ^2 



Al^{vp) ~ 2 dcos^ cose pe(vp cose) ^^^(^^^^^ + • (^50) 
In our energy range Eq. (A7) becomes 



Pe{v)-nUl^v, (B51) 



and so 



AUv,) ^ -i 4 f^v, i |ln(^ ) + 1 ^ . (B52) 
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b. Scattering Part 



The electrons in the hot plasmas that we consider have such large velocities that their 
scattering off the projectiles is quantum mechanical. This is described by Eq. (10.41) of 
BPS which gives 
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With E = ^TTLpVp <^ mpT/mg, the damping constant in the exponent PemeVp/2 is now 
small, not large as it was before. Hence the exponential may simply be replaced by unity, 
and we encounter the integrals 



and 



Hence 
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(B55) 

(B56) 



c. The Sum and a Rough Approximation 

The sum of the dielectric part (B52) and the scattering part (B56) gives 

E:$>T, meE/mj, <T: or T<£;<^T: 

me 



4n \ 2n J 3 



(B57) 



Figure 19 compares this high-energy approximation with the exact result. Figure 20 shows 
that the high and low energy approximations are quite similar. 

Again, to the rough, logarithmic accuracy that produced Eq. (B45) for the ions, we now 
have for the electrons 



An \mp Te 



Note that this electron contribution has the leading temperature dependence given by the 
factor ~ T^^^^ and thus increases as the temperature is lowered. This is in marked 
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FIG. 19: The coefficient Ae (blue curve) compared with the (red curve) high-energy approximation (B57). 
The plasma is equimolar DT with rie = 1.0 x lO^^cm"^ and Tg ~ ~ lOkeV. 

contrast with the corresponding ion contribution given by Eq. (B42) or Eq. (B44) which, 
to within logarithmic accuracy, is independent of the temperature. The ions dominate at 
low projectile speeds as shown in Eq. (Bl), and their contribution at the low speeds also 
behaves as and so also increases as the temperature is lowered. On the other hand, as 

noted immediately below, the electrons greatly dominate at very high projectile speeds with 
a result that is completely independent of plasma temperatures. These remarks provide a 
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FIG. 20: An a particle projectile moving in a equimolar DT plasma with Te — Tj — 10 kcV and Ue = 
1.0 X lO'^^ cni^'^. The (blue curve) low energy approximation (B15) lies above the (black curve) exact result 
while the (red curve) high energy approximation (B57) lies below the (black curve) exact result. Because 
= 2Kg for our equimolar DT plasma, the two approximate forms (B15) and (B57) differ only by a factor 
of two inside the logarithm, and this leads to only slightly different slopes. 
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qualitative description of the stopping power behavior in a plasma. 



4. E ^ nipT/me'- Electronic Contribution 

The high velocity limit in this case has already been calculated by BPS in Eq. (10.43), 
which we simply quote here: 

Vp ^ oo : 

This limit is mostly academic, since the system enters the relativistic regime at these high 
velocities. 



5. Energy Cross Over 

As we have made explicit, the energy loss to the ions in the plasma dominates at low 
projectile energies while the loss is to the electrons at high projectile energies. Here we shall 
estimate the crossover point, the projectile energy at which the two types of loss mechanisms 
arc comparable. Wc shall find that this occurs at a projectile energy that is much greater 
than a typical plasma temperature T, and so we will assume the limit E' ^ T in estimating 
the crossover point. 

For the ions, the E ^ T result (B42) reads 

V y2 I \epeiKe 2 J ' 2} ^ ^ 

This holds provided that 

^-2 _ /47r%y _ / y 2E ^ ^ ^^^^^ 

To put the total ion contribution in a convenient form, we again define a "total ion squared 
plasma frequency" by 

replace the ion charge inside the logarithm by a typical value Cj, and write m^j ~ 77T.p/2 
to approximate the total ion contribution by 
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Note that the only temperature dependence in this result is within the electron Debye 
wave number inside the logarithm. Hence the result only weakly depends upon the plasma 
temperatures. 

A reasonably good approximation for the crossover projectile speed Vp — Vc should be 
obtained by equating the ion result (B63) to the electronic result (B57) which we repeat 
here using 
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In equating the ion and electron approximations (B63) and (B64) we use the crossover energy 
defined by 

^ " (B65) 



to obtain 
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It is important to note that this crossover point only depends upon the electron temperature 
Tg. The ion temperature is of no relevance here. 

Note that the results that we have obtained provide an approximate form for the total 
A coefficient as a function of the energy E — mpVp/2, A{E) — Ae{E) + Ai{E), namely 



A{E) = XE^''' 
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where 
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To return to assess the validity of our approximation for the cross over energy, we examine 
equimolar DT plasmas traversed by alpha particles of mass nip — m^, charge Cp — 2e, and 
initial energy Eq — 3.54 MeV produced by DT fusion. In numerical terms for this case with 
the electron temperature and the crossover energy Eq measured in keV, and the electron 
number density rig measured in cm"^, the crossover relation (B66) appears as 
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FIG. 21: The solution of the crossover condition (B69) as a function of the electron temperature for an 
electron number density = 1.0 x 10^^ cm^"^. The straight line is a fit to these points, with Ec = 51 Tg. 
Similar results obtain for the other densities, with the results presented in Eq. (B70). 

As shown in Fig. 21 for one electron number density, we find that the crossover energies 
for different electron densities Ue are nearly linear functions of Te given by 



C-TeX ( 
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He = 1.0 X 10^^ cm ^ , 

TT-e = 1.0 X 10^^ Cm~^ , 

= 1.0 X 10^6 cm-^ 



(B70) 



Figure 22 shows that for energies below 10 keV or so, these linear relations break down by 
about 10 percent. 
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FIG. 22: This figure displays the lower temperature region of the previous Fig. 21. Here the linear relation 
can deviate from the points which are solutions of the crossover relation (B69) with an error on the order 
of 10%. 

Appendix C: The G Function Simplified 



Here we turn to the definition (4.41) of G{Ti,Te; Eq) in order to reduce it to a more 
manageable form. For convenience, we repeat tliis definition liere^: 



G{Ti, Te] Eq) 



dEF{E) e 



-S(E) 



wfiere 
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We recall that G gives a contribution 

A£;: 



So 



G{Ti,Te\ Eq) . 



Since the prefactor multiplying G involves a temperature difference that is at most only 100 keV and the 
energy Eq is typically 3.5 MeV, this prefactor is less than about 3 %. Hence to within an accuracy of a 
few tenths of %, we need only compute the pure number G to an absolute precision of 0.1 . 
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and 

^^^^-^ {TAm ■ ^^^^ 

First we note that if E' < Eq, the theta function in the curly braces is unity. Hence we 
can make use of the sum rule (4.7) to write 

GiZ, T,; Eo) = ^1(7;, Te; Eo) + G^iZ, T,; Eo) + Gs{T„ T,; Eo) , (C4) 

where 

G.iT.,T,,E.)^l dEEiE)e-si^^l f ^/ dE"V^^e-s^-\ 

(C5) 

roc r^F/ P+S(E') pE' 

(C6) 

and 

roc pE irp/ „+S(E') poo 

(C7) 

First we show that G3 may be neglected. For the very last pair of integrals in G3, since 
the energies E' and E" are larger than £"0 S> , Tj, the electron contribution to the Ah 
functions dominate, and so 

S{E") - S{E') = ^ {E" - E') . (C8) 

J- e 

This is a very large number unless E" is near E' . Hence, with corrections that will be of the 
very small order T^/ Eq, we have 
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and so, again since the electrons dominate the Ah functions in the high-energy regions that 
appear here, 

poo rp pE JTpf 1 

G3(T..T.^o)^-X^ ,E-Ame-^'^>^ l^jj^y (CIO) 
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There is really no need to go any further in the evaluation of G3{Ti, T^] Eq) since it has the 
exponentially small factor exp{— with E > Eq. In this region, as we have noted, the 
electrons dominate and so exp{— ~ exp{— £'/Te}. Even for an electron temperature 
as high as 35 keV and for a DT fusion alpha particle with Eq — 3.54 MeV, this factor is 
exp{-100} ~ 4 X 10-^1 

For the evaluation of G2, it is convenient to define 
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dE"y/E^Jfe-'^^"^ 



To isolate the leading pieces, we shall write 
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and integrate by parts. This will provide an extra explicit factor of a plasma temperature 
T in the numerator, thereby yielding a small quantity. 

The final double integral in the triple integral (C6) defining G2 now appears as 
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where the ellipsis represents the series resulting by further partial integrations. As we shall 
see, the second term in the last line of Eq. (C13) is already negligible, and so are these 
omitted terms. The approximate equalities in Eq. (C13) neglect lower limit terms since they 
result in exponentially small quantities from the remaining integration over E in Eq. (C6) 
because of the factor exp{—S{E)} with E > Eq. Here, to within very good accuracy. 
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since the integral defining H{E) has long since converged to its limiting value at E — Eq. 
Hence, 
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Here, since at large energies the rate of energy variation is of order !/£', 
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{TA{E,)) ^ T_ ^ ^^^g^ 
Eq A{Eq) Eq 



in which T is a typical plasma temperature. The ratio T/Eq is at most a few percent for 
the plasma configurations that we consider, and thus it is a good approximation to replace 
the curly braces in Eq.(C15) by unity. 

Recalling the definition (C3) of F(E) and then using the relation (C12), we obtain 

- OT) C"^ {^^wmp-} -p(-p(^)-s(^o)i> 
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As before, we have the estimate 
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E \ {TA{E)) 

Here again T represents a typical plasma temperature, and since the integration region starts 
at £■ = £'0, we have T/E < T/Eq, Since the factor in the curly braces in the last hne in 
Eq. (CI 8) is just the factor in the curly braces in the first line in Eq. (CI 7), we see that the 
last line in Eq. (C17) is of order T/Eq times the first hne, and thus gives a correction on the 
order of a few percent. We have found that, to within corrections of a few percent, 

G.(T„Te;Eo) - (C19) 

The accuracy of the analytical approximation (C19) for G2 has been confirmed to this 
precision by direct numerical evaluation of its definition (C6). 

In summary, Eq. (C4) expresses the G function in three parts. The first part Gi involves 
a triple integral that must be evaluated by numerical computation. This evaluation is 
simplified because, with the partition that we have made, the regions of integration that 
appear in Gi are restricted to the finite interval < E < Eq. For the second part G2, the 
approximation (C19) is sufficiently accurate for our purposes. The remainder G3 is very 
small and we may simply set 

Gs{Z,T,;Eo)^0. (C20) 
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